Dimensional crossover of spin chains in a transverse staggered field: an NMR study 
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Heisenberg spin- 1/2 chain materials are known to substantially alter their static and dynamic prop- 
erties when experiencing an efTective transverse staggered field originating from the varying local 
environment of the individual spins. We present a temperature-, angular- and field-dependent ^^Si 
NMR study of the model compound BaCu2Si2 07. The experimental data are interpreted in terms 
of the divergent low-temperature transverse susceptibility, predicted by theory for spin chains in co- 
existing longitudinal and transverse staggered fields. Our analysis first employs a finite-temperature 
"Density Matrix Renormalization Group" (DMRG) study of the relevant one-dimensional Hamilto- 
nian. Next we compare our numerical with the presently known analytical results. With an analysis 
based on crystal symmetries we show how the anisotropic contribution to the sample magnetization 
is experimentally accessible even below the ordering temperature, in spite of its competition with the 
coUinear order parameter of the antiferromagnetic phase. The modification of static and dynamic 
properties of the system due to the presence of a local transverse staggered field (LTSF) acting on the 
one-dimensional spin array are argued to cause the unusual spin reorientation transitions observed 
in BaCu2Si207. On the basis of a Ginzburg-Landau type analysis, we discuss aspects of competing 
spin structures in the presence of magnetic order and the enhanced transverse susceptibility. 



PACS numbers: 75.10.Pq, 76.60.-k, 75.40.Cx 

I. INTRODUCTION 

In systems of reduced dimensionality, thermal and 
quantum fluctuations are the source of physical phenom- 
ena without analogues in ordinary 3D materials. Since 
the original studies of Mermin and Wagner on the role 
of thermal fluctuations,— it is known that in the isotropic 
spin-S" Heisenberg model the divergent number of low- 
energy thermal excitations (i.e. spin waves) suppresses 
any long-range ordered phase for dimensions d <2. Due 
to the reduced dimensionality, the magnetic properties 
of spin arrays such as spin chains and ladders are fur- 
ther influenced by quantum effects that are masked in 
common 3D-materials. At T = 0, despite quantum cor- 
rections present in the antiferromagnetic case, magnetic 
order develops in the ground state of the 2D bipartite 
Heisenberg model. Due to the statistical analogy of a 
quantum d-dimcnsional antiferromagnet at zero temper- 
ature with a purely classical magnet at finite temperature 
in d -I- 1 dimensions;^ quantum effects, particularly effec- 
tive in the antiferromagnetic case, imply the absence of 
long-range order in quantum Heisenberg antiferromag- 
nets for d= 1. 

In real bulk materials, however, the same concept of 
one-dimensionality is ill-defined. For instance, it has 
often been shown that the effective dimensionality of 
certain materials can strongly change upon variations 
of temperature^ magnetic field^ or energy scaled be- 
ing probed. As a consequence, the interpretation of 



the experimental observations requires an interpolation 
between limits of different effective degrees of freedom. 
This is particularly true for the specific case of spin- 
1/2 antiferromagnetic (AF) Heisenberg chains realized in 
well-known model materials such as KCuFs,— Sr2Cu03,— 
or copper pyrazine dinitratci^i^ Their high-temperature 
disordered phases, well described by the strictly one- 
dimensional Luttinger Liquid concepti^ are, at low tem- 
peratures, replaced by a magnetically ordered state in- 
duced by weak interchain interactions. The need to inter- 
polate between different dimensionalities suggested the 
introduction of the class of so-called quasi-lD materials. 
The degree of one-dimensionality is usually measured in 
terms of the ratio between the Ncel temperature, marking 
the onset of the 3D magnetic order, and the intrachain 
exchange interactioniiir— Values of this ratio much be- 
low one signify the preservation of the ID character of 
the system. Nevertheless, ID quantum fluctuations con- 
tinue to be effective in the 3D domain too, i.e., below 
the Neel temperature. This is evident from both the 
strongly reduced saturation moment in the ordered state 
as well as from the deviations with respect to excitations 
predicted by the standard spin-wave theory] ^^d^ Obser- 
vations of the evolution of the properties of a physical 
system across the phase diagram, in regions character- 
ized by different effective dimensionalities, has motivated 
a number of theoretical and experimental studies, dedi- 
cated to the phenomenon of dimensional crossover. 

The physics of an assembly of isolated spin- 1/2 quan- 
tum chains upon increasing the interchain couplings is 
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rather well understood i^ The impact of such a per- 
turbation on chains experiencing a local transverse stag- 
gered field (LTSF), however, is still an open questionji^ 
A spin- 1/2 quantum chain in a uniform field and a con- 
comitant LTSF is usually modeled by the Hamiltonian>i^ 

Jf^jY^S,- S,+i - g^i^H J2 S- + (-1) VB^^i ^^ 

i i i 

_ . (1) 

where i is the site-index along the chain direction, J is 
the intrachain exchange coupling constant, /Ub is the Bohr 
magneton, g is the electronic gyromagnetic ratio, z and 
X are the directions parallel and perpendicular to the 
externally applied field H , respectively. The locally in- 
duced transverse staggered field is = cH , with c as a 
constant. 

The interest in spin chains with an LTSF was first 
triggered by experimental studies on copper benzoate 
Cu(C6H5COO)2-3H20^ probing the compound's 
magnetism via susceptibility, nuclear magnetic res- 
onance (NMR) and electron-spin resonance (ESR) 
measurements. In zero magnetic field {H = 0), copper 
benzoate was known to be just another realization of 
a spin-Y2 Heisenberg chain, with an exchange constant 
J/fce ^ 18 K and the onset of a canted antiferromagnetic 
order at ^ 0.8 Kj^ Due to the relatively small intrachain 
exchange coupling, copper benzoate was considered as 
a favorable material for studying the field-dependent 
properties of the ID spin-Y2 quantum Heisenberg model. 
However, the field-dependent bulk properties were soon 
found to be incompatible with corresponding theoretical 
results;^ Incommensurate soft modes at field-dependent 
reciprocal space positions of the excitation spectrum are 
expected for spin-Y2 antiferromagnets in an applied field. 
The first inelastic neutron scattering (INS) experiment, 
intended to verify these theoretical predictions, was 
performed on copper benzoate.— Field-driven incom- 
mensurate modes were indeed found at the expected 
positions in reciprocal space. Due to an unexpected 
field- and orientation-dependent spin-gap A, scaling 
as iJ^/^, these modes were not soft, however This 
surprising INS observation was soon interpreted as the 
result of the particular character of the g-tensors and 
Dzyaloshinskii-Moriya (DM) interactions. Along the 
chain direction, the non-diagonal tensor components 
change sign from one Cu site to the next and the same is 
true for the DM D-vector. Both can be mapped onto an 
effective LTSF and thus to Eq. ([T]). Formally it turns out 
that, in the absence of interchain perturbative terms, the 
effective low-energy theory of the LTSF Hamiltonian is 
given by the quantum sine-Gordon (SG) modelj^i^i The 
SG model is one of the few non-linear problems which 
benefits from an exact solution. Hence, the low-energy 
theory provides analytical expressions for the physical 
quantities in the temperature range A < T <^ J when 
gUsH <C J. 

In the last 15 years the magnetic properties of a large 
number of both organic and inorganic compounds have 



successfully been described by the LTSF model. Specif- 
ically, NMR and ESR studies focused mostly on the 
joint detection of the temperature-dependent longitudi- 
nal (oc {Sf){T)) and transverse (cx {Sf){T)) local mag- 
netization. Examples are the cases of copper pyrimidine 
dinitrate^^'^e and of BaCu2Gc207^. 

Quite generally, any real material, unless directly ex- 
cluded by symmetry arguments, will develop an arbitrar- 
ily small LTSF, captured in the H± term of Eq. ([T|), when 
an external field H is applied. This field-induced LTSF 
is expected to compete with that developing due to inter- 
chain interactions below the Neel temperature (T < Tn). 
The resulting physics in the limit where both effects are 
of similar magnitude is currently not clear.— One of the 
few experimental studies where both perturbations on 
the ID spin-Y2 Heisenberg model were strong enough to 
be experimentally visible is an INS experiment on CuGl2 
• 2(dimethylsulfoxide) (CDC)i^ Spin excitations mea- 
sured at 40 mK showed that the competition leads to 
the opening of a gap at a non zero value of the stag- 
gered field H±, rather than in zero field, as the scaling 
relation A oc would imply. Besides the particular 
case of CDC, all other ESR and NMR studies of materi- 
als modeled by Eq. ([T]) were usually limited to the range 

r>TN. 

By exploiting the high sensitivity of nuclear magnetic 
resonance to local magnetic fields, we present a ^^Si NMR 
study of BaCu2Si207. We demonstrate how the local 
magnetization in a chain, adequately modeled by Eq. ([T]) , 
develops above Tn and how it evolves when T < Tn . The 
latter evolution could be monitored thanks to the com- 
plete decoupling of the local field-induced magnetization 
from the order parameter which is related to the onset 
of a standard, temperature-driven second-order magnetic 
phase transition. 

In section |TT] we briefly summarize relevant results of 
previous research on BaCu2Si207. Our NMR data are 
presented in section IIIIl while a detailed discussion of 
the analysis is reported in section IIVI Finally, section |V] 
offers a summary and the main conclusions of this work. 



II. BaCuaSiaO?: SUMMARY OF PREVIOUS 
RESEARCH 

BaCu2Si207 crystallizes in the orthorhombic space 
group Pnma {D}^) with lattice constants a = 6.862 A, 
6 = 13.178 A and c = 6.897 A^ Each unit cell con- 
tains eight Cu^"'' ions, which determine the material's 
magnetic properties, and eight silicon atoms, each family 
being equivalent by local symmetry. Both types of atoms, 
together with the respective closest-oxygen-atom configu- 
rations, are depicted in Fig.[TJ The copper- ion spins form 
zigzag chains along the crystallographic c axis. Early 
zero-field INS studies of BaCu2Si2 07 at T > Tn revealed 
a gapless spinon continuum. A fit to the lower boundary 
of the spectrum provided an intrachain exchange value 
Ji — 24.1 mcVj^ At the same time, elastic neutron- 
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FIG. 1. a) Schematic view of the crystalhne unit cell of 
BaCu2Si2 07. Red and black numbers identify the copper- 
and silicon atom sites, respectively. The small red spheres 
indicate the oxygen sites. Zigzag Cu'^"'" spin chains run along 
the c axis. 



scattering measurements found evidence of a long-range 
antiferromagnetic order below Tn = 9.2 K, with an or- 
dered Cu^+ moment at saturation of 0.15 ^b, coUinear 
with the c axis.— 1^ These results qualify BaCu2Si207 to 
be among the best realizations of a spin- 1/2 Hciscnberg 
chain system, characterized by a low ordering tempera- 
ture and by a large intrachain exchange coupling. For 
a global overview of a material classification the reader 
is referred to Table 1 in Ref. [ill Over the years, stud- 
ies of BaCu2Si207 proceeded along two main directions, 
both related to the present work: i) the study of the spin- 
reorientation transitions^^— and, ii) the investigation of 
the lD-to-3D crossover)^ 

Regarding the first topic, experiments employing neu- 
tron scattcring j^^i^^ ESR^^ and ultrasound techniqueSf24 
were carried out in externally applied magnetic fields of 
varying strengths and orientations. These studies pro- 
vided evidence of a number of phase transitions related 
to the realignment of spins (and spin-flop transitions for 
an externally applied field H || c) at T < In spite 

of serious efforts, the cause for these transitions is not yet 
clearj^^i^ In the following we show that NMR data pro- 
vide direct evidence of the presence of an enhanced trans- 
verse susceptibility below Tn, offering a simple explana- 
tion for the observed spin reorientations. We suggest that 
non-trivial phase diagrams below the magnetic ordering 
temperature may appear naturally in anisotropic quasi- 
ID antifcrromagnets exhibiting a strong reduction of the 
ordered moment. 

As for the second topic, early studies have shown that 
BaCu2Si207 is also a model spin-chain compound for in- 
vestigating the crossover from quantum-spin ID dynam- 
ics to semi-classical 3D (spin-wave) dynamicsi^ii^ Zero- 
field INS experiments probing the range below the Neel 
temperature revealed how the presence of weak inter- 
chain interactions J2 <Si Ji, produce an effective LTSF 
originating from neighboring chains, which leads to a 
confinement of massless spinon excitations at T < Tn. 
Therefore the situation is not equivalent to that of non 



interacting ID chains in a staggered field. The additional 
J2 interactions induce dispersive excitations with a mo- 
mentum transverse to the chain direction. Their energy 
vanishes at well defined positions in reciprocal space, co- 
inciding with the Bragg peaks of the ordered structure. 
For BaCu2Si207 in an applied magnetic field, intrinsic 
anisotropics lead to a field-induced LTSF causing a spin- 
gap to be present already above Tn. The possibility of 
directly measuring the field-induced gap via INS is slim, 
but important consequences are expected for the local 
static magnetization. Since the latter is accessible via 
NMR, we took advantage of this unique opportunity to 
study the lD-to-3D crossover in a quantum-spin chain 
with an LTSF modeled by Eq. 



III. NMR EXPERIMENTAL RESULTS 

The NMR measurements were carried out on a 4 x 2 x 2 
mm'^ single crystal of BaCu2Si207. Since the a and c lat- 
tice parameters are roughly the same, the correct identi- 
fication of the two crystallographic directions is not triv- 
ial. Field-dependent magnetization measurements on the 
same sample with H || c could detect the two expected 
spin-fiop transitions, hence confirming the correct iden- 
tification of the crystal axes. 

The ^^Si NMR lines were measured as a function of 
temperature at two different magnetic fields, with the ex- 
ternal field being applied along the crystallographic a or 6 
axis. Spectra above the antiferromagnetic transition were 
recorded with a standard spin-echo technique, while at 
lower temperatures the line intensities were established 
by a superposition of frequency sweeps as described in 
Ref. 113. The NMR spectra reported in Fig. [2ji show the 
height-normalized shapes vs. temperature, measured at 
7 T for both crystal orientations. In the paramagnetic 
phase (T > Tn) single lines are observed, which split 
into two in the ordered regime. While the results are 
qualitatively the same for both field directions, quanti- 
tative differences are obvious. With the field along the 
b axis, the NMR line shifts towards higher frequencies 
(by ca. 30 kHz) upon lowering the temperature. A much 
smaller shift is observed when H || a and likewise the line 
splitting for T < Tn is considerably reduced. The signal 
width, instead, is approximately the same in both cases 
and is essentially unaffected accross the entire tempera- 
ture range in the paramagnetic regime. 

The monotonic decrease of the resonance frequency 
with increasing temperature in the paramagnetic phase 
(sec Fig. ^jp) is quite surprising. A spin-Y2 Heisen- 
bcrg chain is expected to display a broad maximum in 
the magnetic susceptibility x{T) at T^ax — 0.64Ji^ 
This maximum, also known as the Bonner-Fisher peak;^ 
should be located at T^ax — 180 K for BaCu2Si207 
(J = 24.1 meV). This was indeed reported in previous 
works^ and confirmed by our low-field magnetization 
data (see Fig. \5!p). Since the NMR line shift and the 
magnetic susceptibility are generally proportional, the 
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FIG. 2. a) Temperature dependence of NMR ^^Si resonances with an applied field H of 7.02 T along the crystallographic b 
and a-axes. b) Relative shift of the peak positions as a function of temperature, c) Direct measure of the order parameter of 
the antiferromagnetic phase, reflected in the difference of the positions of the split lines. The mean values in the ordered phase 
are indicated by black dots. The zero frequency marks the resonance frequency of ^''Si nuclear spins in a standard reference 
sample such as Si(CH3)4.— 



absence of any maximum in Fig. [5}d for T > Tn is a 
striking feature, indicating a fundamental difference be- 
tween the microscopically probed local field (via NMR) 
and that reflected in the macroscopic susceptibility. 

Once the lineshape maxima in the 3D ordered regime 
(T < Tn) were identified (see Fig.[2jD), the differences be- 
tween the peak positions, Av, were evaluated. The cor- 
responding values for the two different orientations are 
shown in Fig. [2J:. We recall that Av is ultimately pro- 
portional to the order parameter of the phase transition 
and hence it can be used to monitor the transition. The 
resulting increase in frequency splitting between 10 and 5 
K can be fitted by a power- law A ■ (Tn — T)^. Although 
the chosen temperature range is too broad to really reflect 
a truly critical regime, the exponent {3 = 0.263 ± 0.004 
is very close to /3 = 0.25, the value obtained from zero- 
field neutron diffraction data.— Figure [2j: also confirms 
that the monotonic increase of the order parameter for 
T < Tn does not depend on field. Indeed, for H || b, prac- 
tically the same ordered moment at saturation is found 
for fiaH = 3.98 T and 7 T. 

The postulated collinear antiferromagnetic order for 
T < Tn is known to have a zero-field saturation moment 
of 0.15 fiB and an easy axis which coincides with the 
c direction.— 1^ As shown in Fig. [21 the line positions 
reflect the NMR response to a magnetic held oriented 
along the a or 6 direction, respectively, i.e., perpendicu- 
lar to the easy axis c. A standard collinear antiferromag- 
net with the held applied perpendicular to the easy axis 
exhibits a constant magnetization in the ordered phase. 
The NMR lines are thus supposed to split symmetrically 



with respect to their common relative shift ^ clearly at 
variance with our observation. In fact, the average posi- 
tions of the two maxima, indicated in Fig.[2lD by full black 
dots, contrary to expectations, are observed to decrease 
with decreasing temperature for T < Tn- Spin- wave cor- 
rections to the constant magnetic susceptibility below Tn 
for antiferromagnets ordered collinearly along a direction 
perpendicular to the applied fleld show at most an in- 
crease of the longitudinal magnetization with decreasing 
temperature. This is due to zero-point spin fluctuations 
affecting the ordered moment.— However, according to 
the data presented in Fig. [SJ:, this correction is modest 
in BaCu2Si207. 



The 29Si NMR resonances of BaCu2Si207 depend 
strongly on sample orientation, a trend which is partic- 
ularly conspicuous in the ordered phase, where the posi- 
tions and shapes of the NMR lines are sensitive to even a 
small degree of misalignment. In order to avoid problems 
with data interpretation due to misorientation, a study 
to establish the orientation dependence of the NMR lines 
was carried out by mounting the sample on a two-axis 
goniometer, suitable for NMR experiments at cryogenic 
temperatures4^ The results for two temperatures above 
Tn, 20 and 200 K, are reported in Fig.[3l Once the a and 
b axes were identified, the sample was rotated such that 
the direction of the externally applied field was kept in 
the ab plane of the crystal lattice. 
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FIG. 3. Angular dependence of the NMR lines at 20 and 200 
K. The orientation of the external field fioH =7.02 T varies 
within the crystalline ab plane. 

IV. DATA ANALYSIS 
A. Origin of the staggered field in BaCu2Si207 

As already mentioned in Sec. HI the zigzag geometry of 
the Cu^+ spin chains in BaCu2Si207 provokes electron 
anisotropies which may strongly affect the physics of the 
chain system in case of an externally applied magnetic 
field. Two dominant contributions to the anisotropy 
originate either in off-diagonal components of the gyro- 
magnetic tensor g (alternating in sign along the chain 
direction), and/or in spin-orbit effects in the Cu-O-Cu 
superexchange path along the chain. The latter is also 
known as the Dzyaloshinskii-Moriya (DM) interaction4i 
To compare the LTSF model of Eq. ([1]) with experimental 
data, an estimate of these two contributions to anisotropy 
has to be made. A gap in the spin-wave excitation spec- 
trum in the ordered phase was observed in zero-field INS 
measurements and attributed to two-ion anisotropy ef- 
fects with an energy scale of ^ 0.4 meV^^ while an ad- 
ditional mode at a lower energy of ^ 0.17 meV was ob- 
served by ESRi^ Based on the local symmetry of the 
intrachain Cu-O-Cu bond, a DM D-vector lying almost 
in the ab plane, with unit vector components [0.86, 0.51, 
0.07] was suggested!^ 

By making use of the crystal symmetry, we apply gen- 
eral space-group operations {Ra,Ta} to a pair ij of cop- 
per sites interacting via oxygen superexchange along the 
c-axis. Rotations Ra and affine transformations are 
related to the symmetry operation a. The original pair 
is transformed into a new set and the local environment 
transforms accordingly. The configuration of the various 
D vectors in the unit cell can be established from the 
transformation rule^i 

where D transforms as an axial vector under the applica- 
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FIG. 4. Upper panel: Pattern of the field-induced LTSF in 
BaCu2Si2 07. Red (blue) lines correspond to the case of a 
longitudinal field applied along the fe(a)-axis. The LTSF ori- 
entation of sites 3,4,7 and 8 is not affected by the change 
of field orientation. Black arrows show the direction of the 
DM vector according to the symmetry rule in Eq. © (see 
text). Lower panel: Local magnetization pattern of the order 
parameter in the magnetically ordered phase at T < Tjv and 
H= 0. The same kind of structure is realized with a moderate 
external field in the ab-plane. 

tion of the rotation Ra- With this rule the full pattern of 
alternating D vectors, depicted as black arrows, halfway 
between the relevant Cu sites in the upper panel of Fig. IH 
can be derived. We note that the a and b components of 
the D vector have alternating signs when moving along 
a given chain in the c direction, or when moving between 
different chains in the a direction.'^ The effect of the 
DM interactions between sites i and j, if not forbidden 
by crystal symmetry, can be taken into account via the 
following spin Hamiltonian4i 

^^DM = D,j-(S, xS,). (3) 

This contribution can be mapped onto a local transverse 
staggered field, H^, via a rotation in spin spacei^ For 
small |D|/J ratios, as is the case for BaCu2Si207, the 
local transverse field at site i can be approximated by: 

ii± ^ ^D,, , X g,rH, (4) 

where gf is the uniform (diagonal) part of the local g- 
tensor valid at site i and H is the externally applied field. 
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The term in Eq. (|4]) represents the second type of the 
LTSF components outhned above. 

In addition, strongly orientation-dependent, high- 
temperature magnetization data were interpreted as in- 
dicating a strong anisotropy of the local g-tensori^ A 
direct measurement of the g-tensor components is in 
general possible via ESR experiment. In the case of 
BaCu2Si207, however, this is hampered by the broad- 
ening and the loss of intensity of the ESR absorption in 
the paramagnetic phase, providing at best an estimation 
of g-factor gt = 2.11 ± 0.07 and = 2.0 ± 0.1,^ The 
evaluation of the local (7-tensor is additionally compli- 
cated by the strong in-chain exchange interaction, lead- 
ing to the exchange (or motional) narrowing of the res- 
onance line.— Consequently, differences in the g-factor 
cannot be resolved as long as the corresponding Zeeman 
splittings are smaller than the exchange energy. Since 
J/fce — 200 K, this is clearly the case here, even in very 
strong fields. As a working hypothesis we assume that 
the (7-factor anisotropy is mostly determined by single- 
ion effects, which reflect the local configuration of oxygen 
atoms. Figure [T] shows that each copper ion is located at 
the center of a tctrahcdrally distorted CUO4 square. Be- 
cause of this distortion, the oxygen-to-oxygen distances 
of opposite O atoms differ by about 2%. If we neglect 
this detail, the local point group of the CUO4 unit is 
D2d- An arbitrary (7-tensor is then invariant under all 
the symmetry operations of the point group D2d- In the 
local reference frame of a CUO4 unit, the tensor adopts a 
uniaxial form g^'^— diag(5i, gi, gg), with two of the prin- 
cipal axes oriented along the two directions at 45 degrees 
from the square's diagonals and the third axis parallel to 
their cross product. The transformation matrix relating 
a Cu04 unit to the crystallographic unit cell may now 
be obtained. Selecting the copper atom 1 in Fig. [1] this 
transformation leads to: 



n.u 

9i 



0.2351 + 0.759 0.31(51-59) 0.33(51-59) 
0.31(51-59) 0.8651 + 0.1459 0.15(59-51) 
0.33(51-59) 0.15(59-51) 0.8451 + 0.1559 

.(5) 

where the subscript i of g'^''^ refers to the site index. Since 
the eight copper sites are equivalent under the allowed 
symmetry operations, we can obtain the tensor 5'^''^ for 
each of them. The matrix in Eq. ([5]) reveals that the 
components 5,^'^ and 5f'^ are roughly equal (in qualita- 
tive agreement with ESR experiments) , reflecting the co- 
inciding high-temperature magnetization tails, measured 
with a field along the b and c direction, respectively. The 
qualitative behavior of the magnetization is also shown 
in Fig. 1 of Ref. |3l|. 

By considering the 5-tensor for each site i in the unit 
cell, and by using Eq. (j4]) for the DM contribution to the 
local field, the expected LTSF pattern for an external 
field applied along the a or 6 axis can be established. This 
pattern is the key for understanding the NMR results. In 
our case, whenever the external field lies in the ab plane, 
the LTSF H]^ is found to be parallel to the c direction (see 
Fig. d]). In order to exploit this favorable configuration, 



the field orientation H || c has not been addressed in the 
present work. With the external field in the ab plane 
and forming an angle 9 with the b-axis, we can calculate 
the ratio H]^ /H. For instance, by considering the copper 
sites 3 and 5 we obtain: 



^1 f 3,2 a , 3,1 • o\ 1 
-Jf = [-91 cos6' + 5i' sm6'j + 

jj (^-Da5i'^ cose* -Db5j'^ sin ( 

H± / 3.2 a I 3,1 ■ , 

-j^ = [9i cos6' + 5i sin( 

^ (Dagl'^ cos 9 - Dtgl'^ sin 9 



(6) 



The configurations of the LTSF, resulting from the ori- 
entation of an external field along the crystalline b (a) 
direction, are indicated by red (blue) arrows in the up- 
per panel of Fig.^l The terms in ([5]) are to be inserted 
in Eq. ([T]). These two patterns are found to be equiv- 
alent to the Lie and L2C irreducible representations of 
the magnetic-structure space, as previously established 
in Ref. l46l They consist of the following linear combina- 
tions: 



Lie — Sic 

L2C — Sic 



S2C 
S2C ■ 



Sac 

S^c 



Sic 
Sic 



S^c 
Ssc 



Sqc 

Sec 



Sjc 
Sjc 



Ssc, 
Shc, 



(7) 



with Sic the component of the local magnetization at the 
i-site along the c-axis. Note that the products LicHb and 
L2cHa are symmetry invariants?^ i.e., they are combina- 
tions of the irreducible representations (IR) of the mag- 
netic structure which transform according to the trivial 
representation (the ID representation consisting of 1 x 1 
matrices containing the entry 1) of the little group of the 
propagation vector fc.' 



B. Static magnetization of a spin-Y2 chain in a 
transverse staggered field 

As mentioned above the local magnetization at a cop- 
per site i, denoted as S*, is given by a uniform and a 
transverse component, such that = SJ^ + S*^, locally 
induced by an external uniform field H„ — 5j"H and a 
staggered field H^. In the following we fix the conven- 
tion that S* has a saturation value of 1/2. In order to 
recall the general results already known for and S^, 
and to present our new results, it is useful to introduce 
the following reduced units: 



h* 
h* 



J ' 



J 



Xl.^ 



d{Sl) 



T 

^ 7' 



(8) 
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FIG. 5. a) Results of finite temperature DMRG calculations 
across an extended range of T* based on Hamiltonian ([1} and 
valid in the small h\ limit, in comparison with known analyti- 
cal and tabulated results (see text for details), b) Comparison 
of the calculated xl a-s obtained analytically or via DMRG 
at T* < 0.4. c) SQUID-magnetometry data measured at 1 
T with the field applied along the a and b crystalline axes, 
respectively. Solid lines are fits using Eq. (|lip . dashed lines 
are predictions for an isotropic uniform chain with no LTSF. 



It has been show n^^i^^ that by using this notation and 
setting both h*^^ ■ and /i* = 0, the susceptibility Xu has a 
peak at T* ~ 0.64 and a zero-temperature limit of I/tt^. 
In Fig. [5^ we reproduce the temperature dependence of 
X*, calculated and tabulated by A. Kliimper in Ref. IH. 
The transverse field h*j_ opens a gap in the excitation 
spectrum which, for /i^ -C 1, scales asi^ 

J = 1.78- (/il)^ • f-lof 



(9) 



For small fields and A/J ^ T* <^ 1 analytic field- 
theoretical results for x*± are availablci^ In the chosen 
reduced units it reads: 



0.2779 log'/=(T 



(10) 



The high h*^ limit has been treated by previous "Density 
Matrix Renormalization Group" (DMRG) calculations.— 
No complete and unbiased numerical result is yet avail- 
able for the susceptibility at small magnetic fields across 
an extended range of T* . In Fig. [5K,b wc fill this gap with 
results of DMRG calculations^"— for chains with 100 
and 200 spin sites, respectively, and compare them with 
the analytical result of Eq. (jlOp . including or omitting 
the logarithmic correction. In our figure, the subscript a 
stands for u or _L, respectively. Without the logarithmic 
term, the main difi'erence between the analytical and the 
numerical result is, as expected, at low temperatures. In 
order to sec how much a transverse field afi'ects the lon- 
gitudinal uniform magnetization, we followed I. Affleck's 
approach^^ and computed numerically the total deriva- 
tive —dF/dhl^ = niu of the free energy of the system 
obtaining: 



(11) 



with the parameters c* as given in Eq. ([6|) . Two separate 
DMRG runs were employed to calculate the uniform and 
staggered susceptibilities independently. Subsequently 
equation (jlip was used to obtain the value of m„. The 
possibility to compute the two quantities (Su) and (S*^) 
separately, allowed us to make use of symmetries to re- 
duce the Hilbert space in the DMRG simulations. 

In order to compare our simulation results with the 
data for BaCu2Si207, the temperature dependence of the 
magnetization f/um„ was measured along the a- and the 
fo-axis of the crystal. The data, measured with a SQUID 
magnetometer, are displayed in Fig. [5]:. Because of the 
small moment, particular care was taken in the choice 
of the sample-holder material that would cause at most 
a small magnetic background signal. With the model 
given by Eq. (jlip we obtain good agreement with the 
experimental data down to approximately 20 K. The de- 
parture of the solid lines from the points is most likely 
due to approaching the onset of magnetic order. As pre- 
vious authors^ we also tried a fit by imposing the value 
c* = (dashed lines). The resulting discrepancies are 
obvious. From the high-temperature tails of m„(T) we 
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extract values between 2.19 and 2 for g^'^ and g^'^, re- 
spectively. The latter differ considerably, but are more 
realistic than the correspon ding values between 2.5 and 
2.2 quoted in Refs. and yd. Useful information can 
be extracted from the fit parameters c'. For the macro- 
scopic uniform magnetization to„, the sign change of 
(sec Eq. ini) is irrelevant. Thus in Fig. [SJ;, and ct are 
the corresponding parameters for fields H along the a- 
or 6-direction, respectively. First of all we consider the 
field- induced and angle-dependent spin-gap in Eq. (0). 
We obtain = 0.61 meV and = 0.78 meV. These 
gaps, when expressed in A/fce units, are both of the order 
of 10 K, i.e. close to the temperature where the magnetic 
order sets in. This explains why an activated behavior 
of the spin-lattice relaxation time has not been observed 
in our previous NMR worki^ A proper estimate for the 
DM parameters Da.b can be obtained by solving Eq. ^ 
with gf^''^ from (O and the fitted Ca,b- We obtain the 
values Da — 0.94 meV and Di, ~ —1.2 meV. Remark- 
ably, the ratio Da/Dt ~ 1.25 is close to 1.68, the value 
obtained from purely geometrical considerations Ex- 
perimentally it turns out that the NMR response at low 
temperatures is dominated by the transverse magnetiza- 
tion {S^), i.e., the diverging susceptibility x*l emerging 
from the DMRG calculation in Fig. [5}d. The study of 
this quantity via NMR and its fate below Tn is the main 
topic of the rest of this paper. 



C. Modeling the NMR lines 

Having established the contributions to the local mag- 
netic field at the Cu sites, we now proceed to study their 
influence on the ^^Si NMR- line data. In our case the 
local magnetization experienced by the silicon nuclei is 
dominated by the externally applied field, and the con- 
tribution due to the sample's magnetization is only of 
second order. For this reason the resonance frequency 
^^Wfc of the ^^Si nucleus k (fc = 1, . . . , 8 — see Fig. [1] for 
the notation) can be written asi^ 



29 



H 

m 



NN 



(l + a,).H + ^n-S* + ^A^S' 



(12) 

where 7 is the ^^Si nuclear gyromagnetic ratio, <7k is the 
orbital shift tensor, is the dipolar tensor which couples 
the silicon nucleus k to the copper atom i, and A], is the 
relevant transferred hyperfine interaction. The dipolar 
sum in Eq. can be calculated directly. This was 
done by fixing the Cu-spin arrangement resulting from 
the LTSF configuration shown in Fig. |4]and by including 
the contributions from the copper atoms within 50 A 
from the considered silicon site. The sum of hyperfine 
interactions runs over the four nearest neighbor (NN) 
copper sites >^ Given and for the silicon nucleus 
k, the relevant tensors for the other silicon sites can be 
obtained by allowed symmetry transformations, meaning 
that tensors referring to the various silicon nuclei are not 



mutually independent. If a symmetry operation of the 
space group {Ra, Tq,}, applied to the silicon site r^, brings 
it to {i?a,Ta}rfc = Tk' (and consequently the copper site 
Vi to r-i'), the hyperfine tensors are given by A^, = ■ 
Aj. ■ Ra- For example: 
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(13) 



The four NN copper atoms surrounding the silicon atom 
located at site k = 1 are i = 3,4,7,8 (see Fig. [ij. 
Contrary to dipolar interactions, the components of the 
transferred-hyperfine tensor are a priori unknown. For 
calculating directly the 0-depcndent component /i^!'^^ of 



the uniform local field at site fc = 1 or 2, parallel to g^H, 
Op, Af = 6^, Al = c„ and A? = d„ and 

m-u and m* = 
(I13L we obtain: 



we define Af = 
we denote m = 
notation of Eq. 



and Af — u.^ 
g^'^{S^). By using the 
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(14) 



where the plus (minus) sign refers to fc = 1 (2). In (fT4|) we 
take into account that the sample's uniform longitudinal 
magnetization and the externally applied field may not 
be collinear due to a possible ^-tensor anisotropy. By 
definition, g"^ = {g^' sinO)^ + (g^' cos^)^ and = o^ -I- 
+ Cfj_ + d^. The reason for picking the Si sites fc = 1, 2 
for describing the relevant NMR lineshapes is evident 
from the contribution of the 
resonance frequency. We get 

( sin 9, cos 



from the contribution of the transverse field h^'^^ to the 
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(15) 



In the paramagnetic phase the following relations always 
hold by symmetry: 



m'^ = — TO^ — nC' ~ —nV^. (16) 

By denoting G3/6 = 03/6 - 63/6 ± C3/6 - (^3/6, we obtain 
for the transverse local field in the paramagnetic phase: 



1/2, 



,3/5 



(G3 sinfl ± Gq cos( 



(17) 



It turns out that by considering any other of the silicon 
sites, the only two possible orthogonal local fields are 
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those given by (|T7l) . For an apphed field along the a or 
the b axis these two local fields coincide. From Eq. © 
it follows that if 6* = 0° (field along 6), = -m^. The 
same is true if = 90° (field along a) (see Fig.U]). A single 
narrow line is thus expected for 9 = 0° and 9 = 90° , while 
two lines are expected in an intermediate angular range 
and at temperatures exceeding Tn- This is indeed the 
case, as already shown in Fig. |31 

We conclude this section with two principal results. 
The first concerns the prediction for the angular-, 
temperature- and field-dependent relative NMR line shift 
Aw in the paramagnetic (T > Tn) regime due to the 
transferrcd-hypcrfine and orbital interactions: 

''Aco,/, =ll^ [Y,gl' + Y,gf + {Y.gf - Y,g\^) cos 20 

±{Y2gf +Yig\^)sm29] 

(G3 sm9±Ge cos 9) 

+ ■^[(71 + (J5 + {(J5 - ai) cos29 ±2a2 sin 26*] | , 

(18) 

where a symmetric orbital-shift tensor a has been in- 
troduced. The consequences of Eq. p8| are discussed 
in the following section. The model behind Eq. is 
ultimately independent of the exact geometry of the hy- 
pcrfinc couplings; the qualitative result does not change 
even if, for instance, the parameters and were zero. 

Next we discuss the situation in the ordered regime, 
below Tn. Here, due to the second-order phase transi- 
tion, the symmetry of the system is spontaneously bro- 
ken. The adopted order reflects one of the irreducible 
representations of the magnetic structure. Its product 
with the applied field is, however, not necessarily an in- 
variant upon symmetry transformations. Consequently, 
in the ordered regime the lines are expected to split, even 
when the field is applied along the main crystallographic 
axes. From previous zero-field diffraction studies^ and 
following the conventions in Rcf. |46| . it is known that the 
representation chosen by the spin system is Tg , coUincar 
with the c axis (see lower panel of Fig. |4]). By calculating 
the local field employing Eq. ()17|) . the NMR lines split 
below Tn according to: 

h'^j^ ^ = mi^{5y^(G3 sm9 ± Gg cos 6*). (19) 

Since in the Lq representation m"^ = is always valid, 
Eq. ()19p suggests that the lines should coincide if 6* = 0. 
This is not the case, however, if all the possible silicon 
sites are considered. In the ordered phase, the local field 
at the sites k = 1,2,4,7 is the same. Also the sites 
k ~ 3,5,6,8 experience the same field, but the latter 
differs slightly from the former. This explains the line 
splitting shown in Fig. [2] for both field orientations with 
respect to the crystal axes a and b. 

The parameters G3/6 = aa/e-^a/e-ca/e + c^a/e 7^ ^3/5 
are not directly accessible by experiment. Nevertheless, 



Eq. (|T9)) offers the possibility to average out the Lg con- 
tribution to the NMR shift below Tn, thus providing a 
direct access to the components m and to* of the local 
magnetization in the ordered regime. 



D. Comparison between theory and experiment 

As just explained at the end of the previous section, the 
average NMR line positions at T < Tn (shown as black 
dots in Fig.[5jD) are independent of the contribution of the 
Lq representation and reflect the influence of the LTSF 
and the uniform magnetization. This holds true even if a 
dipolar term is added to Eq. ([T8)) . since the average NMR 
line position is not affected by the expected symmetrical 
dipolar splitting below Tn. 

With the external field in the ab plane, the LTSF is 
always coUinear with the c axis. For T > Tn the trans- 
verse magnetization induced by the LTSF contributes to 
TO„ in the form of Eq. pT|) . For T < Tn both m and to' 
are still present. Since however m (in Fig. [5j;) is weakly 
temperature dependent, the strong variation of the av- 
erage NMR shift for H || 6 at T < Tn (see Fig. Et) is 
dominated by a contribution related to the Li(T) repre- 
sentation. Thus, below Tn, NMR allows us to reveal the 
effects of the interaction between the representations Li 
and Lq. We return to this issue after considering first the 
T > regime. 

The relative shift captured in Eq. includes the 
anisotropic orbital-shift tensor. Its components along the 
main crystal axes are usually determined via the so-called 
Clogston-Jaccarino plot, where the NMR line shift is 
plotted versus the corresponding susceptibilityj^i^ This 
approach requires a sufficiently broad temperature range, 
in which the NMR shift mimics the sample's magneti- 
zation. Because of the large Cu-0 exchange coupling 
in the BaCu2Si207 chains and the weak ^^Si NMR sig- 
nal at elevated temperatures, a reliable estimate of the 
orbital-shift components was not possible in this way. By 
assuming this tensor to be symmetric and temperature- 
independent, we released the parameters cti, 02 and (T5 
and the hyperfine couplings. In this way the whole 
data set could be fitted with a single set of parameters. 
The quantitative temperature dependences of to and to' 
were established by using the calculations described in 
Sec. IIVBI inserting the values of the LTSF as obtained 
from the fits in Fig. [SJ:. 

In Fig. [nh.-d we display the result of the fits, as well as 
the individual contributions to the local magnetization, 
as a function of the angle 9 and of temperature. In these 
figures the individual contributions to the total shift (red 
curve) of the NMR lines caused by the local longitudinal 
and transverse magnetization and by the orbital shift are 
highlighted as blue, black and green curves, respectively. 
For hqH = 3.98 T, only the global fit is presented. Also 
shown are the temperature independent contributions of 
the orbital shift (broken lines in Fig.[Bl; and d). It may be 
seen that the temperature dependence of the shift due to 
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FIG. 6. Detailed comparison of the microscopic model captured in Eq. [18] with the experimental ^®Si NMR line positions for 
T > Tn- In each panel the green, black and blue lines represent individual contributions to the fit related to the orbital shift, 
the transverse staggered and the uniform longitudinal magnetization, respectively. a,b) Angular dependence of the positions 
of the line maxima with H in the ab- plane, measured at 20 K and 200 K (6 axis corresponds to 6 — 0). c,d) Temperature 
dependence of the NMR shifts measured with the field applied along the b (left) or a (right) axis. The zero frequency marks 
the undisturbed resonance frequency of ^^Si nuclear spins. 



the longitudinal magnetization is weak. The transverse 
component |m^| is small at 200 K but it grows signifi- 
cantly at low temperatures. The data were fitted in the 
range 20 to 230 K. The relative shift of the NMR hnes (at 
very low fields with respect to saturation) scales linearly 
as a function of field. 



The fit parameters we obtain are Yi ~ 0.014 T/fiB, 
Y5 ~ -0.16 T/fiB, Ge = -0.0752 T/fiB and G3 = 0.1286 
T/hb- Since the fit parameters Y2 and Y4 are not lin- 
eary independent, we could fit only their combination 
123^'^ + Y^g^'^ = 0.53 T//xb. The computed dipolar ten- 
sor components, to be inserted in Eq. are of the 
order of 0.02-0.04 T/fiB- The orbital shift values dis- 
played in Fig. [S]are of the order of 150 ppm. 



E. Competing spin structures 

Employing the same classification of representations as 
introduced in Ref. |H, the antiferromagnetically-ordered 
phase in zero magnetic field is related to the Lq repre- 
sentation and the corresponding order parameter. In the 
previous section we provided evidence for an enhanced 
transverse magnetic susceptibility even in the ordered 
regime. This enhancement is characteristic of quasi-lD 
chains in an LTSF and we argue that it is the reason 
for the unusual spin-reorientation transitions that are 
observed in BaCu2Si207. The microscopic approach re- 
quires considering the effects of the lD-to-3D dimensional 
crossover in specific features of the magnetic properties. 
In case of chains with no LTSF, this was done with a com- 
bined mean-field and "Random-Phase Approximation" 
approach)^ 

Here we tackle the problem with a Ginzburg-Landau 
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(GL) expansion^ of the free energy ^ close to Tn- Al- 
though this phenomenological approach neglects fluc- 
tuation effects, it has the advantage of retaining the 
exchange-energy contributions to the susceptibility of the 
ordered phase. The mean-field approach also includes 
interactions between different, possibly coexisting, order 
parameters. 

We start by constructing symmetry invariants of the 
little group of the k-vectoiM In BaCu2Si207, even if ex- 
posed to an applied field, a commensurate antiferromag- 
netic structure with fc = is realized, leading to a little 
point group which coincides with D2h- We call the ^- 
component of L^^ , with the /3-IR; clearly the product 
13 ■ fjL = 'iN (with = 8 as the number of equivalent cop- 
per sites in the unit cell) . In terms of a GL free-energy ex- 
pansion over all possible order parameters, a phase tran- 
sition will occur whenever one of the coefficients Af} oi the 
quadratic term ApL'jj^ changes sign. Since the igc repre- 
sentation is the one realized in the magnetically ordered 
regime in zero field, we can write that Aq = e^iT — T^) 
{e% > 0). All the IRs L^g^ of the magnetic structures in 
BaCu2Si207 are one-dimensional4^ It is therefore easy 
to construct invariant combinations of the Lp^ since the 
representations of the powers of these terms, which have 
to transform according to the trivial representation, re- 
main one dimensional. The expansion in Eq. (|20p is based 
on the physics discussed in the previous sections of this 
paper. Retained are the terms containing Lg^, related to 
the zero-field magnetic order, representing the LTSF, 
and H the external magnetic field. We will limit our 
considerations here to the case of an LTSF pattern Li e 
which is realized for H || 6. A more complete analysis 
will be published separatelyi^ The relevant expansion in 
powers of Lg/J: and reads: 

= 00 + AeLl + AiLj + BqLI + SigL^L? 
+ B'.eiU ■ + D{U ■ Uf + D'n^Ll + aall^ 



H2 

Stt' 



(20) 



with (j)Q as the value of the free energy in the paramag- 
netic phase in zero field. 

The Ai coefficient is positive above the Neel temper- 
ature, refiecting the absence of a spontaneous symmetry 
breaking related with the Li order parameter. On the 
other hand we assume Ai to be small in the vicinity of 
Tn since, as indicated in Fig. SI Lie differs from the low- 
est energy-state configuration Lgc only by the mutual 
orientation of the spins in the neighboring chains, which 
are relatively weakly coupled. The temperature depen- 
dence of Ai is assumed to be linear in the vicinity of Tn: 

The fourth order term for Lgc fixes the magnitude of 
the main order parameter below the transition; as re- 
quired, Bq > 0. The terms with prefactors Big and 
B[q arc crucial in our discussion, because they describe 



the exchange competition of the field-induced order Li 
and the spontaneous order Lg. Microscopically, these 
terms arise from the simple idea that both the main or- 
der parameter and the induced order parameter involve 
the same local spins, eventually along the same crystal- 
lographic direction. The Bie coefficient is expected to be 
positive in order to enhance the energy cost for the coex- 
istence of these two magnetic structures. Finally the term 
related to B[q defines the preferred mutual orientation of 
the two order parameters by means of the scalar product 
between the them. From the expansion in Eq. ([20)) alone 
it is not possible to predict whether a collinear (Sjg < 0) 
or a transverse {B[q > 0) spin configuration is realized. 

The terms related to D and D' describe interactions 
between the longitudinal magnetization and Lg- The 
terms with prefactors describe the orientation of the 
zero-field order parameter. Reported results of neu- 
tron diffraction^^ and antifcrromagnetic resonance^ im- 
ply that Oa > Oh > 0. The term ab is responsible for 
the fact that the magnetic structure Lie is induced by an 
external field applied along the b axis. The powers of H 
completing the expansion in Eq. (|20p provide a full de- 
scription of the effects of the g-tensor anisotropy on the 
longitudinal magnetization. This may be seen by recall- 
ing that: 



-— = M. 

(9H in 



(21) 



Minimizing over the components of Li for H || 6 and 
assuming a zero- field collinear antiferromagnetism Lgcj 
we get Lia ~ Lib = and: 



Lie 



-oibHb 



2Ai{T) 



1 



^Ll 



(22) 



With a similar reasoning we obtain the longitudinal mag- 
netization Mb along the fe-axis, in the former notation 
denoted as m„ : 



Mb = -2D(tl-Le)L6c 



2D'LlHb 



o-bLic + ixp + lb)Hb- 
(23) 



The last two equations deserve some discussion. 

Equation (|22p captures the temperature dependence of 
the transverse staggered magnetization. Microscopically, 
the increase of Lie upon cooling above Tn (Lgc = 0) 
is due to the divergent transverse susceptibility of a ID 
spin-Y2 quantum Heisenberg chain in an LTSF. This situ- 
ation is modeled by the decrease of Ai{T) on cooling (i.e. 
with e'l''' > 0). It also predicts a decrease of Li upon the 
growth of Lg at T < Tn, as observed in the NMR data. 
Microscopically, the change of regime upon decreasing 
temperature, from a divergent transverse susceptibility 
(characteristic of a spin chain) to a progressive compe- 
tition between the field-induced magnetization pattern 
and the zero-field order parameter, is argued to be a di- 
rect consequence of the dimensional crossover from ID to 
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3D of a chain in an LTSF. The clear experimental iden- 
tification of how ID physics affects the static magnetiza- 
tion properties even below Tn is the new result emerging 
from the present study. Below we address the question 
of how these anomalous properties for T < Tn can ex- 
plain certain spin reorientation transitions observed in 
BaCuaSizOy. 

By analyzing Eq. ([23|) we note that the first term is 
zero for an easy axis (c axis in our case) orthogonal to 
the direction of the applied field (along the b direction). 
The second term, instead, provides corrections to the 
constant magnetization predicted by the standard mean- 
field theory below Tn- Microscopically it can be related 
to a semi-classical contribution of spin-waves4^ From the 
magnetization data in Fig. [5}: it may be concluded that 
D' < 0. Next we single out a constant paramagnetic con- 
tribution, with prefactors Xp E^^id 7b, which is related to 
the magnetization of the ideal spin- 1/2 Heisenberg chain 
at the Neel temperature. We note that also in this case 
a contribution from the staggered transverse susceptibil- 
ity affects the longitudinal magnetization data, M{T). 
While Lie scales as abHt,, the contribution to Mf, scales 
as alHb, in full qualitative agreement with the micro- 
scopic approach. In magnetization measurements along 
either the a or the b axis, the contribution of the stag- 
gered magnetization matters, but it is not as outstanding 
as in NMR measurements, where both Alf, and Ld are 
revealed. 

The enhanced transverse susceptibility accounts very 
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FIG. 7. Data representing relative ^''Si NMR line shifts in 
the vicinity of the ordering temperature, compared with the 
classical predictions of Eqs. (|22[1 and (|23|l . Although fully 
classical, the GL type approach grasps the competition be- 
tween the spin structures described by Lie and Lq. The hy- 
perfine parameters and the orbital shift are taken from the 
analysis for T > Tn (see text). The average NMR frequency 
is shown below the transition temperature. Black, blue and 
green line show single components of the fits due to Lie, Mb 
and a, respectively. 



well for the observed spin reorientations. For example, 
with the applied field along the 6-axis, such a transition 
occurs at Hsi ~ 7.8 Ti^ It is caused by a sudden change 
of the easy axis from the c- to the a-direction. Below Tn 
the quasi one-dimensionality extends itself in the form 
of a field-induced transverse susceptibility. It favours a 
field-induced spin alignment Lie which competes with 
the zero-field order parameter Lqc- The higher the field, 
the larger is the energy cost to sustain this arrangement 
[captured by the term Bie of Eq. Substituting 
Eq. dm into (^0]) yields an expression which depends on 
Lg only. A spin reorientation is then expected as the 
result of competition of the field- dependent anisotropic 
corrections with the conventional anisotropy of the order 
parameter at a field: 



Hsr — 



2M 

at 



B' 



(24) 
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It can be shown that, for B[q > 0, Eq. (PU)) also accounts 
for two spin-reorientation transitions when H || c. The 
inclusion of the staggered field pattern described by the 
representation L2c [see Eq. ^] could similarly account 
for the phase transition at H || Based on formu- 
las (|22p and (|23p . we now attempt a comparison with 
the experimental data in the temperature range T < 20 
K, where the microscopic ID model does not properly 
describe our results. Using the transferred hyperfine pa- 
rameters determined in section llVDl we compare the ^^Si 
NMR line shift monitored for a field Hs,- = 7.02 T ori- 
ented along the 6-axis with the GL approach, postulating 
Lec = [(Tn -T)/Tn]^ and /? = 0.5 (we arbitrarily set the 
zero-temperature limit of Tgc to unity). This is shown 
in Fig. [T] To obtain tentative estimates of the GL-model 
parameter, we first fitted Mi, of Fig. [5J: in the vicinity of 
Tn to Eq. (|23|) and obtained the parameters Xp + 7& ~ 
7.05 X lO"'^ cmu/mol Cu, D' = -7.6 x 10"^ emu/mol 
Cu and the ratio ag/(2^J°^) « 7.01 x lO^^* cmu/mol Cu. 
Notice that the present numbers contain already a pref- 
actor 0.15 in the numerator of Eq. (P^ . corresponding 
to the experimentally reported zero-temperature limit 
of Lqc expressed in Bohr magnetons. Considering the 
decrease of Lie with increasing temperature above Tn 
we obtain e\°^ — 0.051 K~^. Next, with the fixed 
al I (2j4^°'' ) ratio wc could fit the relative NMR peak po- 
sitions, as shown in Fig. [71 and hence determine the pa- 
rameters J^i^ and Sie + i?5.g. The fit shown in Fig. [7| 
was obtained with J^^^ w 3.15 x 10^ emu/mol Cu and 
Sie + B'l^ ss 7.63 x 10'' emu/mol Cu, which correspond 



to a value Lic(Tn) 
7.02 T. 



0.15 ■ ai,Ebl{2Af') 0.035 Ate at 



V. SUMMARY AND CONCLUSIONS 

A detailed analysis of ^^Si NMR data obtained by 
probing single-crystalline BaCuaSiaOr revealed the in- 
fiuence of ID physics into the regime of 3D magnetic 
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order at temperatures below 10 K. In this way the prob- 
lem of weakly interacting nearest neighbor chains, de- 
scribed in the non-interacting limit by the model in Eq. 

could be addressed. Based on a classical Ginzburg- 
Landau analysis it is shown that in this type of com- 
pounds complicated {H, T) magnetic phase diagrams 
emerge. They are caused by the interaction of the trans- 
verse staggered local magnetization, originating from 
magnetic anisotropics in spin-i/2 Heisenberg chains, with 
the effective magnetic field due to the weakly ordered 
spin moments on neighboring chains. We argue that 
the previously established spin-reorientation transitions 
in BaCu2Si207 reflect this situation and can, therefore. 



be understood in this framework. 
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